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ABSTRACT 

We investigate a scenario of photons scattering by electrons within a relativistic out- 
flow. The outflow is composed of discrete shells with different speeds. One shell emits 
radiation for a short duration. Some of this radiation is scattered by the shell(s) be- 
hind. We calculate in a simple two-shell model the observed scattered flux density 
as a function of the observed primary flux density, the normalized arrival time delay 
between the two emission components, the Lorentz factor ratio of the two shells and 
the scattering shell's optical depth. Thomson scattering in a cold shell and inverse 
Compton scattering in a hot shell are both considered. The results of our calculations 
are applied to the Gamma-Ray Bursts and the afterglows. We find that the scattered 
flux from a cold slower shell is small and likely to be detected only for those bursts 
with very weak afterglows. A hot scattering shell could give rise to a scattered emis- 
sion as bright as the X-ray shallow decay component detected in many bursts, on a 
condition that the isotropically equivalent total energy carried by the hot electrons is 
large, ~ 10 52-56 erg. The scattered emission from a faster shell could appear as a late 
short 7-ray/MeV flash or become part of the prompt emission depending on the delay 
of the ejection of the shell. 
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1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are a cosmological phenomenon 
with a huge energy release, fast variabilities and very com- 
plex multi- wavelength light curves. A relativistic outflow is 
unavoidable in order to explain the fast variability and so 
called "compactness problem" (see Piran 2005 for a review). 
According to the standard "Fireball" model, the outflow 
from the GRB central engine has a finite duration and can 
have a wide range in its velocities, thus can be modeled by 
being made of discrete relativistic shells. These shells are re- 
sponsible for the observed 7-rays (via internal shocks) and 
for the afterglow emissions (via external shocks) (cf. Piran 
2005). In this picture, if one shell emits 7-rays, some frac- 
tion of that emission should be scattered by shells behind, 
and the scattered emission would arrive at the observer at 
a different time, with a different flux and possibly at a dif- 
ferent photon frequency. Detection of the scattered photons 
would help us explore the properties of the GRB ejecta and 
the late outflow. 

In the present paper we consider a simple scenario, 
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where only two consecutive shells are present: one shell radi- 
ates and the other receives some of this radiation and scat- 
ters it. The two shells can have different speeds and the shell 
that receives and scatters the emission can have an arbitrar- 
ily large time delay in its ejection from the central source, 
but it has to be behind the emitting shell. An observer de- 
tects the primary emission from the first shell and then the 
scattered emission from the second one at a later time be- 
cause of the light-travel time. We calculate the ratio between 
these two emissions' fluxes, the time delay in their arrival, 
the ratio between their frequencies and the ratio between 
their durations. 

Early GRB X-ray observations by Swift have shown a 
"canonical" behavior that presents a puzzling shallow decay 
typically lasting for a few hours (e.g. Nousek et al. 2006). 
This decay phase is poorly understood (see Zhang 2007 for 
a review of current possible models). We will explore the 
possibility that this shallower decay could be due to the 
scattered emission. 

The scattering of the GRB prompt emission photons by 
electrons or dust grains in a dense circum-burst environment 
has been investigated before (e.g., Esin & Blandford 2000; 
Madau et al. 2000; Shao & Dai 2007; Heng et al. 2007). The 
scattering process we consider in this work happens within 
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the GRB outflows, which is a natural outcome of the outflow 
when it has a finite duration and a variable speed. 

The paper is structured as follows. We first derive a gen- 
eral formula for the observed flux from a relativistic shell in 
§2. In §3 we construct the flux and geometrical relations for a 
two-shell model. The formulae for the scattering process are 
developed in §4. Then we elaborate the primary and scat- 
tered emission relations such as time delay, frequency ratio 
and time duration ratio in §5. We present the main result 
- the ratio between the scattered and the primary fluxes - 
in §6. The application to the GRB shallower decay data is 
presented in §7. A faster scattering shell case is discussed in 
§8. We also discuss X-ray dim bursts and X-ray-dark short 
bursts, for which the scattered emission might be easier to 
detect, in §9. Finally, the summary and conclusions are given 
in 810. 



2 EMISSION FROM A RELATIVISTIC SHELL 

Consider a spherical shell moving relativistically with 
Lorentz factor (LF) Y (when the shell is beamed with an 
opening ang le ^ r _1 , it still can be considered as being 
spherical). The surface brightness in the rest frame of the 
shell is e' v , (erg s~ cm~ Hz ~~ 1 sr _1 ), the luminosity distance 
between the observer and the shell is Dl and the radius of 
the shell is R, both of these distances measured in the labora- 
tory frame. The flux density received at a frequency v by the 
observer ahead of the shell, /„, can be obtained by calculat- 
ing the specific luminosity of the relativistic shell. The lumi- 
nosity of the shell is given by f v {4mD\) = e'„, (47ri? 2 )r(27r). 
The last expression includes a factor of T, to take into ac- 
count the "boost" that the photons experience; and a factor 
of 2n, since we assume that the photons are being emitted 
isotropically from the shell, in the rest frame of the shell, and 
we are only interested in the ones reaching the observer. The 
two expressions yield 



/„ = 27re'^r 



R 
D L 



(1) 



3 TWO SHELLS SCENARIO 

Consider now two thin shells being ejected with an half open- 
ing angle of 9j from the central engine. Shell 1 is ejected first 
with LF Ti and, after a delay St, measured in the laboratory 
frame, shell 2 is ejected with LF T2 ■ We assume that shell 1 is 
emitting photons isotropically in its co-moving frame and is 
characterized by an angular-independent surface brightness, 
e' v , , on both sides of the shell. In the laboratory frame most 
of these photons will appear to move in the same direction 
as shell 1 and reach a distant observer, and a few will move 
in the opposite direction, encountering shell 2 on their way. 
These photons will be scattered by shell 2 and then reach 
the observer. See Figure 1 for an illustration. We will use 
primed quantities to specify the co-moving frame where the 
quantity is being measured: unprimed corresponds to the 
laboratory frame, primed (') to the co-moving frame of shell 
1, and double primed (") to the co-moving frame of shell 2. 
Before proceeding with the detailed calculations about 



the scattered emission, we provide simple scaling relation- 
ships between the observed scattered emission and the ob- 
served direct emission from shell 1 by only considering the 
line-of-sight region. 

Let us assume v is the photon frequency of the direct 
emission from shell 1. In the shell 1 co- moving frame, the 
emitted photon frequency is v' ~ due to the relativistic 
Doppler effect (for simplicity we neglect the factor of 2). As 
seen by shell 2 the photon has a frequency of v" ~ u' j- 2 - . If 
shell 2 is cold, the scattering does not change the photon's 
energy. The observed frequency of the scattered emission is 
v s ~ v"Yi. Thus the observed frequency ratio between the 
two emissions is — ~ (-p 2 -) 2 . 

The scattered emission will be observed at a later time, 
because the scattered photon has traveled an extra distance. 
This extra distance is equal to Ri(l— ^j+fcSt, where R\ is 
the distance of shell 1 from the central engine when a photon 
was emitted from shell 1 toward shell 2, and /3i and /?2 are 
shell 1 and shell 2 velocities, respectively (the first part of 
the expression is due to the difference in the speeds of the 
shells, and the second part is due to the ejection delay of 
shell 2). The observed delay of the scattered emission is this 
separation divided by the speed of light « i^k^ ~ T(^) 2 for 
small St, where T is the observed time of the shell 1 direct 
emission since the central engine explosion. 

If L is the luminosity observed directly from shell 1, the 

shell 1 co-moving frame luminosity would be L' ~ -ki ■ In the 

1 1 

shell 2 co-moving frame, shell 1 has a luminosity of L" . Using 
the Lorentz invariance of I^/v 3 we find that the luminosity 
ratio, ^jt, is the frequency ratio to the fourth power, or 

l£r ^ (t^) 4 ^ (£f) 4 - In the shell 2 co-moving frame, r e L" 
is the luminosity of the scattered emission, where r e (< 1) 
is the shell 2 electron's optical depth. Then the observed 
luminosity of the scattered emission is L s ~ T e L"T\. Thus 
we obtain that ~ r e (^-) 6 , which shows that the observed 
luminosity ratio is strongly dependent on the LF ratio. Since 
~ (f^) 2 , this indicates that the observed specific flux 

ratio is ^~r e (£a-) 4 . 



3.1 Incident flux on shell 2 

Now we work on the calculation of the scattered flux on 
more detail. In order to determine the scattered flux from 
shell 2 we will first calculate the incident flux from shell 1 at 
the point of shell 2 that intersects the line of sight between 
the central engine and the observer - we will call this point 
P. To do this, we will use the Lorentz invariance of Iv/v 3 , 
where I„ (erg s _1 cm -2 Hz -1 sr _1 ) is the specific intensity 
and v is the frequency of the photon. 

Let us consider a bundle of rays being emitted at an 
arbitrary point - Q - on shell 1 and directed to the point 
P (see Figure 1). The angle between the line of sight and 
the line connecting the central engine and Q is 6. The angle 
between the line of sight and the line connecting P and Q is 
a. From the Lorentz invariance of I„/v 3 we can obtain the 
relation 



T" 



r 2 (l +/3 2 cosa) 



(2) 



where I„ is the specific intensity of the bundle of rays 
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Figure 1. A simple two-shell scenario geometry. The diagram shows two photons emitted at the same time from point Q on shell 1. One 
photon travels to the observer ahead of shell 1, while the other travels back to point P on shell 2. Note that the figure shows shell f and 
the photons at time t and shell 2 at time t + I, where / is the light travel time from Q to P. 



measured by an observer standing still in the laboratory 
frame and at the position of point P (note that I u is NOT 
the specific intensity of the emission detected by a distant 
laboratory- frame observer sitting in front of shell 1); /"„ is 
the specific intensity of the bundle of rays measured in the 
shell 2 co-moving frame at point P. We can relate I v with 
the specific intensity of the bundle of rays measured in the 
shell 1 co-moving frame at point Q, I' v , , as follows 

' 1 ^ (3) 



F, ITi[l + /3icos(a-0)] 



We need to obtain a relation between the co-moving 
specific intensity, I' v , , and the surface brightness, e' u , , of shell 
1, both quantities in shell 1 co-moving frame. This relation 
is given by 



t, = 



COS T) 



(4) 



where r)' is the angle measured in the shell 1 co-moving 
frame between the photon's direction and the normal to the 
emitting surface (facing shell 2) . r\ can be determined using 
the aberration of light formula (Rybicki & Lightman 1979): 

, cos(a-0)+/3i 
COSV = l +/ 3 lC os(a-0) (5) 

where all the quantities have been defined previously. 

Finally, using formulas (2) - (5) one can obtain the inci- 
dent flux from shell 1 on point P in shell 2 co-moving frame, 
/"// , which is given by 



/"" — j I'v" 



cos a" dQ," — 2ir 



cos a sin a da , 



rll I ± 2 

Jv" = Ke,,! ( jr- 



(1 + /? 2 cos a) 3 d sin 2 a' 1 



+ cos(a - 0)][1 + /?i cos(cv - <9)] 2 



(6) 



where the integral runs from to a" - the half opening angle 
of shell 1 as seen from an observer on point P in shell 2 (the 
subscript "j" denotes the edge of shell 1). 



3.2 Light path geometry 

In order to solve this last integral, we need to understand 
the relation between the angles and how they transform in 
different frames. First, we need to setup the geometry of the 
problem. The photon that is emitted from point Q at time 
t travels a distance I and reaches point P on shell 2 at time 
t + 1 (here, and throughout the paper, we use units in which 
the speed of light is 1). The radius at which the photon was 
emitted from point Q on shell 1 is Ri(t) and reaches point 
P on shell 2 at radius R2(t + l), where R(t) can be obtained 
by R(t) — /3t (all these quantities are measured in the lab- 
oratory frame). It is important to remember that there is 
a time delay in the ejection of shell 2, which will be taken 
into account when using R2(t). The geometry describing the 
light path of the photon gives: 



Icosa = R t (t) - R 2 (t + l) 

I sin a = 7?i(i)sin#, 

from which it can be shown that 
ft 

tana = 



1 - 



f3 2 , hit 



fee 



(7) 
(8) 

(9) 



where we use Ri to denote Ri (t) and we have also used the 
fact that 6 is small, so that sinf? ~ 9. And again, using the 
aberration of light formula we have 



tana — 



r 2 (cosa + /32)' 



(10) 



Also from (9) and (10) we can get a close form expres- 
sion for a": 



tana ' 



(ii) 



Equations (10) and (11) allow us to carry out the inte- 
gral in equation (6). Careful analysis (see Appendix) shows 
that the integrand is weakly dependent on a" and it reduces 
to a constant of order unity. This simplifies /"// greatly and 
we obtain: 
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f"„ = ne' v > (jr^J J dsh>2 a " = ne 'v' f J sir<2 a 'i ^ 

where a" can be obtained from equation (11) by setting 
6 = 6,. 



Recall that before the detailed calculation we estimated 
the observed specific flux ratio under the line-of-sight ap- 
proximation. Eq. (14) is consistent to what we estimated 
earlier, except that the previously ignored shell solid angle 
term is fully considered here. 



4 SCATTERING FROM SHELL 2 
4.1 Scattered flux 

Assuming that shell 2 is "cold" and knowing the flux from 
shell 1 at point P, we can calculate the scattered flux from 
shell 2. We will assume that the photons from shell 1 un- 
dergo Thomson scattering with the electrons on shell 2 and 
that the scattering is isotropic in the co-moving frame. The 
surface brightness of shell 2 in its co-moving frame will be 
given by 

1 



4tt 



where E e is the electron surface density and gt is the Thom- 
son scattering cross section. To obtain the scattered flux at 
the scattered frequency that reaches an observer at a lumi- 
nosity distance Dl one can use equation (1) 



27re„//r 2 



R 2 (ts 



Dr 



and obtain 

fl,s = TT"T2 



R2 it aC at 



ft 



(13) 



where we have r e = E e <TT, the optical depth for electrons to 
Thompson scattering. We also have used Ri{t sca t) to specify 
that the radius of the second shell needs to be calculated at a 
later time, t 3ca t, when the scattering occurfl Using equation 
(12) we obtain 

r =-7 

which allows us to find a ratio between the scattered flux 
from shell 2, /^ s , and the direct flux from shell 1, /„, 





Rl it scat) 


m 





f 

T 



= — T e Sin Of, 



Ri{t s 



Ri(t) 



(14) 



where we have used equation (1) to immediately get f v . It 
is important to note that these two fluxes arrive at different 
times. This will be further explained in detail in the next 
section. Also, the result for r e > 1 is the same as that for 
r e = 1. 



4.2 Shell radii 

All quantities of the flux ratio are known, except the ratio of 
the two radii. To obtain it, we will consider a photon emitted 
on shell 1 that travels along the line of sight to shell 2. The 
light travel time of this photon is given by equation (7) but 
using a = 0, so that 

I = Ri(t) - R 2 (t + l). 

Solving for I we obtain 

= fli(t)-fi 2 (t) 
I + P2 

The ratio of the radii is given by 
Rl{tscat) _ fait - St + 1) 



(15) 



Ri(t) 



Pit 



and using equation (15) and the fact that both shells move 
close to the speed of light, we find 



Rljtscat) _ fh _ fhSt 

Ri(t) ~ Pi 2Ri ' 



(16) 



With this last equation and equation (11), the flux ratio 
given by equation (14) - is fully determined. 



4.3 Time dependence of scattered emission 

For simplicity, let us assume that the emission from shell 1 
is constant and time independent: that it is a box function 
with some finite duration. The time dependence of the scat- 
tered emission will be given by: (i) the time evolution of the 
electrons' optical depth, r e , (ii) the opening angle of shell 1 
as seen by a co-moving observer on shell 2, a'j , and (iii) the 
radius of shell 2, R2(t sca t)- 

The time dependence of R^itscat) and a'j is weak, pro- 
vided that the ejection time delay between shells, St, is small 
compared to Ri. This is the case we are interested in, since 
if St ~ iii, then St would be on the order of hours or days. 
This scenario would invoke a very long lasting activity of the 
central engine, a scenario that we don't want to address in 
this paper. The only time dependence of the scattered emis- 
sion will be given by the time evolution of r e , which goes as 

T e OC R- 2 OC T~ 2 . 



1 We approximate the incident flux on point P to the incident 
flux on any other point on shell 2. The validity of this approxi- 
mation depends on the magnitude of a'.'. If a" is small, the an- 
gular size of shell 2 as seen by an observer co-moving with shell 

2 at the point where the line of sight intersects with shell 1 - let 

us call it A" - is also small, because —i? = " R2 i ta /:? t ^ < 1; then 

this approximation should be good. For our interested parame- 
ter space - determined from the GRB data we are going to use 
- a" is between 0.1 and (see §7.1.2), not very small. Thus the 
approximation overestimates the scattered flux slightly. 



5 PRIMARY AND SCATTERED EMISSION 
RELATIONS 

5.1 Time delay 

Let us assume that two photons are emitted from shell 1 
at the same time. Photon 1 travels directly to the observer 
located ahead of shell 1 and arrives at time T v (p stands for 
primary) . Photon 2 travels in the opposite direction, scatters 
from shell 2 and then travels back to the same observer 
arriving at a later time T a (where s stands for scattered). 
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What is the time delay, T 3 — T p , between the arrival of these 
two photons? 

If we are only interested, as a simplification, in the pho- 
tons along the line of sight, then this time delay will be given 
by equation (15). We only need to multiply this expression 
by 2, to obtain the full time it takes for the photon to travel 
to shell 2 and then to travel back to shell 1. Then, the time 
delay is 

Ri(t)-R 2 {t) 



Ts Try — 2 - 



Further simplification yields 



T v = Rj 



Wi 1 " hl)t ' 



(17) 



5.2 Ratio of frequencies 

In this section we will determine the relation between the 
frequency of a photon emitted from shell 1, v, and a photon 
emitted from shell 1 and then scattered by shell 2, v s (both 
quantities measured in the laboratory frame). For this, we 
will only consider the photons traveling along the direct line 
of sight between the central engine and the observer. 

In the shell 1 co-moving frame, a photon is emitted 
from shell 1 with frequency v' . In the laboratory frame, this 
frequency is measured as 

v = i/ri(i + /?i cose'), 

where 0' is the angle between the photon's direction and the 
x'-axis direction of the shell 1 co-moving frame. Note that 
the a;'-axis of frame K' is always directed to the moving 
direction of frame K' relative to another inertial frame (K 
or K"). Since the photon is emitted to the observer ahead 
of shell 1, 0' = 0, and the previous equation becomes 



i/ = i/ri(i + /?i). 



(18) 



Let us now consider a photon emitted in shell 1 and 
directed to shell 2. The frequency of this photon in the shell 
1 co-moving frame is v' , and in the shell 2 co- moving frame 
is given by v" . They are related by 

V = V V x (1 + /?! COS ), 

where F" (/?") is the LF (velocity) of the shell 1 as mea- 
sured in the shell 2 co-moving frame. These quantities can 
be expressed in terms of quantities in the laboratory frame 

as 

|/3i-&| 



r' 1 ' = rir 2 (i-/?i/3 2 



01 = 



1 - P1P2 ' 



If we assume that Ti > Y2 (or Ti < r 2 ), in the co-moving 
frame of shell 2 it will seem that shell 1 is moving away from 
shell 2 (moving towards shell 2), so that 0" = tt (0" = 0). 
Using the last 3 equations, we can determine 



i/Tir 2 (l-/3i)(l + /32 



(19) 



which holds for both assumptions. The photon will be scat- 
tered by shell 2 by Thomson scattering (there will be no fre- 
quency change in the co-moving frame of shell 2) and then 
will travel towards the observer. The scattered frequency 
can be obtained with 

Vs = i/T 2 (l +/3 2 cos0"). 



Setting 0" w 0, since the photon moves towards the ob- 
server, yields 



v s = i/T 2 (l + /? 2 ) 



(20) 



Finally, using equations (18)-(20) the ratio of the scat- 
tered frequency to the primary frequency is 

2 

(21) 



> 2 / x 2 

r 2 \ ( 1 + & 



ri 



l + /3i 



using the fact that both shells travel close to the speed of 
light. This is consistent with our earlier simple estimation. 
Therefore, a slower (faster) shell 2 will lower (raise) the fre- 
quency of the primary photons. 

5.3 Ratio of observed durations 

Let us assume that shell 1 emits for a finite duration of time, 
At (in the lab frame). An observer located in front of the 
shell will detect that the radiation from shell 1 lasts for AT P , 
given by AT P = At(l — /?i). The radiation from shell 1 will 
also travel back to shell 2 and will get scattered, giving a 
scattered radiation duration of AT S in the observer frame. 
If the first photon from shell 1 is emitted at time t and the 
last one at t+At (in the lab frame) , then we can use the time 
delay equation (17) to find the time delay between the first 
primary photon and the first scattered photon, and the time 
delay between the last primary photon and the last scattered 
photon, respectively. Subtracting these two expressions, we 
find that 

& ~ (ll 2 
r 2 



AT S 



AT V 



l-/3i 



(22) 



which means that the observed duration of the scattered 
emission will be stretched (shortened) by a factor of (jj^j 
for a slower (faster) shell 2. 



6 RESULTS 

6.1 Ratio between scattered and primary fluxes 

We want to write equation (14), the ratio of the fluxes, in 
such a way that we can easily use the available observations 
to test our theory. 

We can use equation (17) to solve for St, the time de- 
lay between the ejection of the two shells. This expression 
is then substituted into equation (11) to get sin 2 a" and 
into equation (16) to get the radius ratio. Finally, we obtain 
the flux ratio in terms of the time delay divided by -S^ as 
follows: 



MT P ) 



(Ts 



2r 2 (9 3 R 1 /2T\ 



02 



1 (T s 



1 

2 + Wx 4ff R 1 /2Tj 



2 / N. 4 

rr) ' <23) 



We then plot this ratio of two flux densities as a func- 
tion of the ^--normalized time delay (the choice of this 

normalization will become evident on the next sections) for 
various values of r 2 /Fi in Figure 2. For plotting purposes 
we assume Ti = 100 and 9j =0.1. Since the exact value for 
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Figure 2. The ratio of the flux of the scattered emission to the flux of the primary emission (measured at their respective arrival times 



and photon energies) as a function of the observed time delay between these two emission components in units of 



Hi 



assuming shell 1 



LF Ti 
ratio. 



100 and shell opening angle 0j = 0.1. The unknown shell 2 optical depth r e is put on the y-axis as an "unit" of the flux density 



r e is unknown, we choose to scale r e into the y-axis of the 
figure, as an "unit" of the flux density ratio. 

Looking at figure 2, we can observe that the theoretical 
flux ratio curves have two regions (this is very noticeable 
when T2 < Ti). In region I the flux ratio is flat, and in 
region II the flux ratio is proportional to the square of the 
inverse of the normalized time delay. If we inspect equation 
(23), we can separate its two regions by: 



•r T s —T p ^ 2r i 



2D j 



«i/2r= 



if 



(Region I), 
(Region II). 



These two regions will be used in our applications section. 
For now, they just provide a simpler theoretical model. No- 
tice that the T2 > li curves are dominated by region II, 
while the Y2 < Ti curves have a combination of both re- 
gions. For the latter, the maximum scattered flux is given 
by region I. 

These results show that the scattered flux from a slower 
shell is small, falls at a lower energy than the primary pho- 
ton's energy and its total duration is larger than that of the 
primary emission. If we have a faster shell, then the scat- 
tered flux from it could be either larger or smaller than the 
primary flux, depending on its ejection time delay, St. If St is 
larger (smaller) than the total observed duration of the pri- 
mary emission, then the scattered emission would appear at 
late times (would be part of the primary emission). In any 
case, the energy of the scattered photons would be larger 



than that of the primary ones, and the total duration of the 
scattered emission would be smaller than that of the pri- 
mary emission, so that the scattered emission would appear 
as a short bright flash. 



6.2 Hot shell 2 

In the last sections we assumed Thomson scattering, but 
we should also look at the possibility that shell 2 could be 
hot, so the scattering mechanism at work would be inverse 
Compton. Shell 2 electrons could be hot in the following 
scenario. Consider that shell 2 is ejected after shell 1 and 
undergoes particle heating, by either internal shocks or mag- 
netic dissipation, at a radius smaller than the radius where 
shell 1 produces 7-ray photons. Shell 2 would experience 
adiabatic expansion and would cool, but the electrons could 
still be hot by the time that the 7-ray photons from shell 1 
reach them. The shell 2 electrons might also cool via radi- 
ation, but the chances still exist that the radiative cooling 
is very inefficient, for instance, when the radiation mech- 
anism is synchrotron-self-inverse-Compton instead of pure 
synchrotron for the same observed typical photon energy, so 
that the electron cooling time could be comparable to the 
delay between the electron heating and the scattering. 

The main difference in the formulas previously derived 
will be in the ratio of the frequencies of the primary and the 
scattered emission. Equation (21) is modified to include the 
inverse Compton boost to the photon energy: 



(24) 



where 7 e is the electrons' thermal Lorentz factor. The the- 
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oretical flux density ratio previously derived will not be 
changed in this new scenario. 

Let us now calculate the isotropically equivalent total 
energy in the shell 2 hot electrons. For this, we need to cal- 
culate the isotropically equivalent total number of electrons 
in the shell from their optical depth 

N e = AtyRI — . 

OT 

In the present scattering scenario, we have Ri fa Ri (we'll 
prove this on section §7.2.2) . Ri can be estimated from the 7- 
ray variability time scale: (explained in detail on §7.1.1). 

Then the total energy in the hot electrons of shell 2 is given 
by 

E e = N e m e c 2 Y 2l , = 16tt f^A T 5 lle m e c 4 ( (25) 

\1 1 J ) OT 

where c is the speed of light. 

7 APPLICATION TO THE SHALLOW DECAY 
COMPONENT IN GRB EARLY X-RAY 
AFTERGLOWS 

In this section we compare the 7-ray burst and X-ray after- 
glow data with the results of the last section to find out if 
the shallow decay in X-ray "canonical" afterglow light curves 
can be due to the scattered emission from the shell(s) follow- 
ing the 7-ray shell. We first consider simply the Thomson 
scattering mechanism, then we turn to consider the inverse 
Compton scattering where the electrons in the scattering 
shell have highly relativistic thermal energy. 

7.1 Thomson scattering 

7.1.1 Data set 

From a Swift GRB early X-ray afterglow catalog presented 
in O'Brien et al (2006), we choose a sample of 10 bursts all 
of which show clearly a "canonical" behavior that includes 
a shallow decay component. We apply our simple scenario 
to these bursts, and assume that: (1) the last 7-ray photon 
that was emitted from shell 1 traveled directly to the ob- 
server (primary emission); (2) at the same time and from 
the same site on shell 1, another photon traveled to shell 2, 
was scattered, and eventually became the last X-ray photon 
of the shallow decay phase (scattered emission). Therefore, 
we will use the ratio of the flux density at X-ray energy at 
the end of the shallow decay and the flux density at 7-ray 
energy at the end of the 7-ray emission. Theoretically, this 
ratio should fit our equation (14) if the shallow decay were 
to have its origin in the scattered emission scenario. 

For the 7-ray photons, the catalog in O'Brien et al. 
(2006) only gives a mean BAT flux. We use this mean flux 
as an approximation to the flux at the end of 7-ray emission; 
the flux density at a specific photon energy is obtained using 
the BAT spectral index Pbat (/„ oc u~ 0BAT ). For the X-ray 
photons, we use the XRT flux at the end of the shallow decay 
and the X-ray spectra index fix [jv oc v~^ x ) to obtain the 
flux density at a specific photon energy. 

In our theory R\ is an unknown parameter, but we 
should be able to extract it from the available data: 



is the 7-ray burst duration Tgo for FREDs (fast rise and ex- 
ponential decay) - those bursts whose light curves are made 
of one smooth or two overlapped pulses; and for bursts with 

multiple spikes in light curves, < Tgo. Note that 

1 IL 1 

is equal to the curvature time scale - the delay between 

two photon's arrival times, one emitted from shell l's vis- 
ible edge and the other from the center of shell l's visible 
region. Thus, the 7-ray variability time scale, if we assume 
it is mainly determined by the curvature time scale, would 

be a good approximation for Therefore, we look up the 

1 n 

7-ray light curves from the Swift archive!; for FREDs, 
is simply Tgo, but for those spiky bursts, we use the dura- 
tion of the last pulse. The data obtained for our sample of 
10 bursts has been organized in Table 1. 

7.1.2 Comparison between observations and theory 
predictions 

For the X-ray shallow decay r2/Fi = 1 is a limiting case, 
since for r2/Ti > 1: (i) the scattered emission from the 7- 
rays would fall at a higher energy, not in the X-rays, accord- 
ing to equation (21), and (ii) would have a smaller duration 
than the 7-rays duration according to equation (22), which 
is not what it is observed - the shallow decay in X-ray light 
curves typically extends up to 10 4 s. 

In Figure 3, we plot our 10 GRBs sample data and the 
results of our theoretical calculations for two cases. For the 
first case, the 10 data points use the XRT flux at 1 keV 
(scattered emission) and the BAT flux at 100 keV (primary 
emission), which corresponds to r2/Ti =0.1 according to 
equation (21). For the time delay, we use (T s — T p ) ^ t en d, 
where t en d is the end time of the shallow decay - we do this 

because Tgo >C t e „d- For we use the method described 
21 1 

in the last subsection. These values are also listed in Table 
1. 

In the same figure, for the second case, the 10 data 
points use the XRT flux at 1 keV (scattered emission) and 
the BAT flux at 10 keV (primary emission). This corre- 
sponds to T2/ri = 1/vTO « 0.3, according to equation 
(21). We consider this case to see if the shallow decay might 
be produced by the scattering of the low energy tail of the 
7-ray emission. Since in this case shell 2 is faster than in the 
first case described above, the time delay in the ejection of 
the shells will be larger. 

The normalized time delay for the sample has a range 
of 2 x 10 2 < J - T3s < 2 x 10 3 . For fiducial parameter values 

V\= 100, T 2 = 10 and S = 0.1, using equations (11) and (17), 
this range has a constraint on a": 0.1 < tana" < 1. Since 
this angle is not very small, the approximation of the ho- 
mogeneous incident flux on shell 2 that was used in deriving 
the scattered flux formula in §3.1 will slightly overestimate 
the observed scattered flux (see discussion in footnote 1). 

7.1.3 Results 

Figure 3 shows that the observed flux ratios of the sample 
are generally (10 3 — lO 4 )^ 1 times larger than the theo- 
retical expectation. The discrepancy would be a factor of 

2 http:/ /heasarc. gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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Table 1. Our sample of 10 GRBs which show clearly an X-ray shallower decay component and the relevant data. Basic data are from 
O'Brien et al. (2006). 



GRB 


050315 


050319 


050713A 


050713B 


050714B 


050803 


050814 


050819 


050822 


050915B 


7-! n 

FbatZ. 

(10 _8 erg cm -2 s _1 ) 


3.2 


0.5 


3.8 


3.2 


1.1 


2.4 


1.2 


0.9 


2.5 


8.6 


q b 
PBAT _ 


1.2 


1 


0.(5 


0.5 


2 


0.5 


1 


1.6 


1.5 


1 


^=100fcev(^90)^ 

/in nT-iA 
(1U /ijyj 


4 


0.8 


12 


6 


0.(5 


4.8 


2 


1 


3.2 


1(5 


(10 uJv") 


63 


8 


48 


19 


(50 


15 


20 


40 


101 


160 


00 


96 


150 


130 


130 


50 


90 


111 


3(5 


105 


10 


w 


30 


40 


20 


130 


50 


90 


144 


3(5 


30 


40 


t j 9 

''end ^ m 

(10 SJ 


1 


3.2 


1 


4 


5 


2 


> 6 


2 


1.3 


5 


z A Ki Vena / 

(10 -11 erg cm -2 s _1 ) 


0.8 


0.8 


1.2 


1 


0.5 


1 


< 0.05 


0.04 


0.7 


0.09 




1.5 


2 


1.3 


0.7 


4.5 


0.7 


1.1 


1.2 


1.(5 


1.5 




1.2 


1.1 


2.1 


1 


0.16 


1 


< 0.06 


0.04 


1.2 


0.1 


tend fc 

«i/(2rjf) — 


330 


800 


500 


308 


1000 


222 


> 417 


556 


433 


1250 


F 1/ =ifceV(tend) 
^i>=100feeV t^bo) 


0.03 


0.137 


0.02 


0.017 


0.026 


0.021 


0.003 


0.004 


0.037 


0.0006 


iWfceV(*end) 
F u=\0keV ( T 9o) 


0.0019 


0.0138 


0.005 


0.0052 


0.0003 


0.0066 


< 0.0003 


0.0001 


0.0012 


0.00006 



a Mean BAT flux. 
b BAT spectral index. 

c BAT flux density at 100 keV at the end of the 7-ray emission. 
" BAT flux density at 10 keV at the end of the 7-ray emission. 
e Duration of the 7-ray emission. 

f Simply T90 for those bursts with one smooth or two overlapped pulses; for those "spiky" bursts, i.e., those with multiple, separated 
pulses, we use the duration of the last pulse. 
9 Ending time of the X-ray shallow decay. 
^ XRT flux at the end of the shallow decay. 
* XRT spectral index. 

XRT flux density at 1 keV at the end of the shallow decay. 
™ Observed time delay between the primary and the scattered emissions in units of See §7.1.2. 



~ 10 4 — 10 5 for a modest value of r e ~ 0.1. This indicates 
that the emission of the shallow decay is too luminous to 
be interpreted simply by the scattering within the two shell 
scenario. 

The same figure shows a smaller discrepancy between 
the sample data and the theoretical curve for a case where 
shell 2 is slightly faster but still not exceeding the 7-ray 
shell's speed. Using a modest value of T e ~ 0.1 would make 
this discrepancy to be ~ 10 2 . 



7.2 Inverse Compton scattering in shell 2 

7.2. 1 Comparison between observations and theory 
predictions 

To compare the fluxes of the X-ray shallow decay with the 
theoretical expectation of the inverse Compton scattering 
scenario, we fix the scattered frequency, v a = 1 keV. With 
equation (24) we can determine the frequency of the primary 
emission. The flux ratio data points in Figure 3 have to be 
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Figure 3. Similar to Figure 1 but with the GRB sample data added. The stars (squares) are our GRB sample data using the ratio 
between the XRT flux density at the end of the shallow X-ray decay at 1 keV and the mean BAT flux density at 100 keV (10 keV). The 
dashed (dotted) line is the result of our theoretical calculation. 



changed because the frequency ratio is changed. The flux in 
the 7-ray band follows /„ oc v^Pbat . With this, we can now 
determine that the data points in Figure 2 will be multiplied 
by a factor of [^-(j^-) 2 7e]~' 3i5AT to account for the inverse 
Compton scattering effect, where v 1 and vx are the 7-ray 
and X-ray photon frequencies, respectively, at which the flux 
densities are used in the data points of Figure 3. 

Since 7 e is an unknown quantity we cannot determine 
where these new points would lie on a plot analogous to 
Figure 3. For this reason, we ask: what is the value of 7 e 
necessary to lower all data points to the theoretically ex- 
pected curve? To answer, let us fix the value of the Lorentz 
Factor ratio, Y%IY\ = 0.1 - this is a reasonable value accord- 
ing to equation (22), since the 7-ray duration usually lasts 
~ 10 2 s and the shallow decay extends to ~ 10 4 s. 

The required values for 7 e are presented in Table 2. We 
present two different values for each burst. We use the BAT 
spectral index (subscript c) and an average between the BAT 
and XRT spectral indices (subscript e), respectively. We do 
this because we are extrapolating the 7-ray flux to energies 
below the observed BAT band and it is unknown if the BAT 
spectrum will behave as a single power law in this region. 



7.2.2 Electrons' energy in shell 2 

We calculate E e for each burst in our sample using equation 
(25). The results are presented in Table 2, adopting Ti = 
10 2 , — 0.1 and r e = 0.1. Note that we have corrected 
for the cosmological time dilation effect. Two values 

21 1 

of energies from the two values of j e obtained in the last 
subsection are given in the table. 

Note that from equations (16) and (17) one can see that 



W = h + TS? - • In our data sam P le 

^ s . ^ 10 3 . Therefore, the approximation made in the 

derivation of equation (25), that R2 ~ Ri for Ti ~ 10 2 , is 
valid. 

7.2.3 Results 

Table 2 shows that the isotropically equivalent total energy 
carried by the electrons of a hot shell 2 is large, ~ iq52-56 
erg. If we take into account the cooling of the electrons via 
adiabatic expansion and/or radiation, the initial total en- 
ergy when the electrons were just accelerated is even bigger. 
The optical depth, r e , would certainty decrease the total 
energy but only by a small fraction. The prompt emission 
from those electrons in shell 2 would arrive at about the 
same time as the scattered emission, and the two emission 
components would have similar durations because both du- 
rations are oc 1/F 2 ,. Depending on the ejecta properties, e.g., 
the ratio of the shell 2 energy to the shell 1 energy, the shell 
2 prompt emission might dominate the scattered emission 
in the light curve. If that is true, the shell 2 prompt emis- 
sion might be a possible origin of the late X-ray flares in 
bursts for which both the flares and the shallower decay are 
present. 



8 FASTER SHELL 2 

In the previous section we assumed a slower shell 2. What 
if shell 2 is faster than shell 1? Based on the formulae we 
have, if T2 > Fi , the scattered emission from 7-rays would 
fall in higher energies, e.g., ~ MeV, not in the X-rays, and 
have a shorter duration than the 7-rays. 



10 Shen, Barniol Duran & Kumar 



Table 2. Calculated values for the required 7 e and total electron energy of shell 2 to obtain the theoretical expected flux ratio. The 
calculations were done with two different spectral indices: the BAT spectral index and the average between the BAT spectral index and 
the XRT spectral index. 







Hi b 


7e* 




^1 


E e J_ 


GRB 


Redshift^ 














to 


[Te"*] 


(erg) 


[r."*] 


(erg) 


050315 


1.95 


10 


34 


5 xlO 52 


23 


3 xlO 52 


050319 


3.24 


9.4 


316 


5 xlO 53 


16 


5 xlO 52 


050713A 




5.3 


1364 


1.4 xlO 54 


95 


5 xlO 52 


050713B 




34 


2404 


1.5 xlO 56 


657 


3 xlO 55 


050714B 




13 


13 


2 xlO 52 


5 


7 xlO 51 


050803 




24 


1931 


6 xlO 55 


547 


1 xlO 55 


050814 


5.3 


23 


25 


2 xlO 53 


22 


2 xlO 53 


050819 




9.5 


10 


1 xlO 52 


13 


1.5 xlO 52 


050822 




8 


21 


1.6 xlO 52 


19 


1.4 xlO 52 


050915B 




10.5 


33 


6 xlO 52 


17 


3 xlO 52 



a References for known redshifts: GRB 050315: Kelson & Berger (2005); GRB 050319: Fynbo et al. (2005); GRB 050814: Jakobsson et 
al. (2006). For bursts without measured redshift, we use the mean redshift z= 2.8 of the Swift GRB redshift distribution (Jakobsson et 
al. 2006). 

^ Cosmological time dilation corrected curvature variability time scale; it is equal to the i^<z in Table 1 divided by (1 + z). 

i 

c Required shell 2 electron thermal LF using Pbat as the spectral index. 

^ Isotropically equivalent total energy of electrons with 7 e calculated in the last previous column for r e = 0.1. 
e Required shell 2 electron thermal LF using ( Pbat+Px ^ as spectral index. 

f Isotropically equivalent total energy of electrons with j e calculated in the last previous column for T e = 0.1. 



It is shown in Figure 2 that, in the F2 > Fi cases, the 
scattered emission is very bright though it decreases with 
increasing time delay. If the scattered emission spectrum 
mimics the power law form of the primary emission spec- 
trum, at some time delay significantly larger than the dura- 
tion of the primary emission, the lower-energy-extrapolated 
scattered flux is close to or even brighter than the prompt 
7-rays. According to Figure 2, for ~ 3 (corresponding 
to ~ 10 according to equation (21)) and a selected ob- 
served time delay (T s — T„)/(^-) ~ 100 (which corresponds 

21 1 

to an ejection delay St ~ 10 3 s), ff, s (T a ) w 0.1[r e ]/„(T p ). Ex- 
trapolating the scattered flux density from v s to v, we have 
fl{T s ) = (^fBAT x .l[r e ]/,(T p ) « fu{T p )\Te] for /W « 
1. For a smaller observed time delay, the flux is even greater. 
That means we should have seen a lagged very short 7-ray 
flash at ~ 10 2 — 10 3 seconds after the burst, provided that 
T e « 1. This case cannot happen because the observations 
have never showed this feature. Even though the flux of the 
scattered emission is smaller for a larger observed time de- 
lay, in order for the delayed scattered 7-ray flash to indeed 
happen but below the BAT flux limit, shell 2 must have an 
extremely large ejection time delay St ^ a few x 10 3 seconds 
(cf. equation 17) which is very difficult to explain in terms 
of the central engine activity. For r e < 1, the F2 > Ti case 
could have happened but the flux would be too small to be 
detected. 

One possible case of T2 > Ii is that the shell 2 ejection 
delay St is small and shell 2 has moved very close to shell 1 
when the scattering happens so that the observed time de- 
lay of the scattered emission is comparable or smaller than 
the duration of the primary emission. The scattered emis- 
sion in the BAT energy range will become part of the ob- 



served prompt 7-ray emissions in time. Future high energy 
(^ MeV) observations (e.g., by GLAST) during the burst 
may be able to determine the existence of this case. 



9 X-RAY DIM OR DARK GRBS 

In §7, our calculations show that the scattered flux is much 
lower than the observed X-ray shallow decay flux. If the 
scattering indeed happens, then the resultant emission must 
have been buried in the shallow decay. Thus, in order for the 
scattered emission to be detected, not only the shallow X- 
ray decay component must be absent, but also the normal 
external forward shock component must be extremely weak 
or absent. In rare cases Swift does observe X-ray afterglows 
without a shallow decay and without the normal forward 
shock decay {ax ~ 1 — 1.3), which we call X-ray dim GRBs. 
They show instead a very steep flux decay {ax 3) and are 
thought to be located in extremely low density regions. The 
steep decay component can be explained by the large an- 
gle emission (Kumar & Panaitescu 2000). GRB 050421 and 
GRB 051210 are two examples (Godet et al. 2006; La Parola 
et al. 2006). Note that GRB 051210 is of the short burst class 
which, according to its compact binary progenitor model, is 
more probable to occur in a low density environment. At 
late times, however, when the large angle emission is low 
enough, neither burst shows any sign of re-brightening atop 
the steep power-law decay. 

Another case of interest is the existence of X-ray dark 
GRBs: short bursts (GRB 050906 and GRB 050925) that 
show no X-ray afterglow detection, only upper limits at ~ 
100 s (Nakar 2007). It is also believed that these bursts took 
place in extremely low density environments. 
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X-ray dim and dark GRBs provide a great opportunity 
to put into work the theory presented in this paper. If we 
assume that there is a late ejecta behind the 7-ray producing 
source, then we can try to constrain its LF by looking at 
these bursts. Since their afterglow emission is so weak (or not 
present), we can ask the question: what are the constraints 
on the LF of the late ejecta, so that the scattered emission 
is present in X-ray dim and dark GRBs, but doesn't exceed 
the actual flux observations or upper limits? We will devote 
this section to answer this question. 

We'll start by using the data from Swifts BAT and XRT 
observations as follows: 



fUTs) < f v X f T (T a 



MT P ) f? AT (T 90 ) 



= 0.01 



fl00kev( T 9°) 



where we have selected v = 100 keV, but we have made no 
assumption on the value of u s only that it should obey equa- 
tion (21). We have also assumed fix = 1, which is consistent 
with Swift observations. 

With the previous inequality and using the two pre- 
viously defined regions of the flux ratio (§6), we can find 
constraints on the LF ratio. In region I, we find: 



and in region II: 



9 J ,_ 1 r 1 , 2 ,1.26/^ 6 T- 1/6 



(26) 



20 



<rf* 



0.1/i/ 4 



T 8 T-n 



Ri/2T\ 



1/2 



l/2 r -l/2 -1/4 



-1 1 1 



(27) 



where we have used f x 



f BAT 



and the convention 



Qx = Q/10 x has been adopted. 

After presenting these last two conditions, we should 
return to the physical picture. The initial assumption we 
made is that the scattered emission should be present in- 
dependently of the burst's data and its region. This means 
that the LF of the late ejecta could be very small, so that its 
contribution to the flux would be also minuscule. Therefore, 
we should discard the lower limit on equation (27). Now, we 
want to constrain the LF of the late ejecta from above; we 
are interested in knowing what's its maximum. The maxi- 
mum value between the two upper limits of equations (26) 
and (27) will be the best and more conservative value to 
choose: 



20 



Hl/21^ 



-iri 



ri 



max < 



o.i/# 4 



1.26/ X7 

T s-Tp 

fli/21^ 



1/6 -1/6 



1/2 



1/2 -1/2 
7 j,-l 1 1,2 7 



-1/4 



Now we are ready to consider the 4 bursts mentioned 
at the beginning of this section and obtain the constraints 
for the LF of the late ejecta. 

Using the last X-ray detection of GRB 051210 (La 
Parola et al. 2006) and the upper limit of GRB 050421 
(Godet et al. 2006), both at ~ 10 3 s, we can obtain the ratio 
between the XRT flux density at 1 keV (scattered emission) 
and the BAT flux density (Chincarini et al. 2007) at 100 



keV (primary emission), f xl . Following the procedure out- 
lined in §7.1.1, the values for J" are also obtainecjf]. We 

K1/21 1 

can do exactly the same for the X-ray dark short GRBs, 
using the upper limits of GRB 050906 and GRB 050925 at 
« 100 s (Pagani et al. 2005; Beardmore et al. 2005; Nakar 
2007). With f xl and the normalized time delay, we find the 
following constraints: 



r, \l/2 1/4 

-1^1,2) V < 




for GRB 050906 and 050925, 
for GRB 050421, 
for GRB 051210. 



The upper limits for the X-ray dark short GRBs provide 
the best constraints, since the upper limits in their X-ray 
flux are very strict. These results imply that, for r e — 1, 
Ti — 100 and 6j =0.1, the late ejecta is very slow T2 < 15, 
but it could be faster if r e decreases. 



10 SUMMARY AND CONCLUSIONS 

We have investigated a scenario of photons scattering by 
electrons within a relativistic outflow. The outflow is com- 
posed of discrete shells. One front shell emits radiation, ob- 
served as the GRB's prompt 7-ray photons. Some fraction of 
the radiation is incident backwards to the shell(s) behind, 
and is scattered isotropically in the local rest frame. The 
scattered emission arrives at the observer at a later time, 
T s , and at a different photon energy, v s , that are determined 
by the LF ratio of the two shells and the time delay of the 
ejection of the second shell. We calculated the flux density 
ratio, i. e., the flux density of the delayed scattered emission 
to that of the front shell's primary emission, as a function of 
the normalized arrival time delay and the assumed LF ratio. 

The calculated flux density ratio are compared with ob- 
servations using a sample of Swift GRB X-ray afterglows 
which show a distinct shallower decay component in their 
light curves, with the motivation to see if the scattering sce- 
nario could be the origin of the shallower decay. The results 
are negative. For Thomson scattering, the flux density of the 
scattered emission is about 10 3-4 T e -1 times lower than that 
of the shallower decay component, where r e is the scattering 
shell's electron optical depth. 

We also consider the inverse Compton scattering sce- 
nario in which the electrons in the scattering shell is hot. 
We find that, in order for the scattered emission flux to be 
bright enough to match the shallower decay component, the 
isotropic equivalent of the total energy carried by the hot 
electrons is large, ~ ]^g52-56 er g rp ne p rom pt emission from 
the scattering shell appears at the same time as the scattered 
emission and with a similar duration. 

In the cases where shell 2 is faster than shell 1, when 
extrapolated to the BAT energy band, the scattered flux 
can be as bright as the emission from shell 1. The delay of 
the scattered emission is determined by the ejection delay 
of shell 2. When the ejection delay of shell 2 is much larger 

3 When following §7.1.1 to determine the normalized time delay, 
there is some uncertainty determining the width of the pulses 
since the data shows statistical noise. For GRB 050421, two pos- 
sible values for the normalized time delay were obtained hence 
two constraints for the LF ratio were derived. The value reported 
here is the more conservative one. 
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than the duration of the primary emission, the scattered 
emission would appear as a late short 7-ray/MeV flash. For 
a small ejection delay of shell 2, the scattered emission would 
become part of the observed prompt emission. The fact that 
no late short 7-ray/MeV flash is detected does not support 
the existence of a late faster shell. 

Lastly, we study the possibility of detection of the scat- 
tering emission in two X-ray dim GRBs - that only show a 
very steep flux decay and do not show either a X-ray shal- 
low decay nor the normal forward shock component - and 
in two X-ray dark short GRBs - that show no X-ray af- 
terglow detection at ~ 100 s. Assuming that there is slower 
moving ejecta material behind the fast 7-ray producing shell 
in these bursts, we find upper limits for the Lorentz factor 
of the late slower material. More sensitive observations of 
X-ray dark short GRBs could provide stronger constraints 
on the presence and properties of slower moving material 
accompanying the fast 7-ray jet in GRBs. 

Almost simultaneous to the appearance of this paper, 
Panaitescu (2007) presents a similar work on the scatter- 
ing of the GRB early emission photons by a late outflow. 
Both papers present the same physics and the main formu- 
lae are consistent. The main differences are that (i) we have 
considered the scattering within the relativistic outflow and 
the photons to be scattered are the prompt 7-ray photons, 
whereas Panaitescu (2007) considers the scattering of the af- 
terglow forward shock photons by a late relativistic outflow, 
and (ii) Panaitescu (2007) has considered a faster second 
shell to be able to explain features like the "shallow decay" 
and the X-ray flares, whereas our focus is mostly on a slower 
second shell. 
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APPENDIX: APPROXIMATION TO THE 
INTEGRAND FUNCTION IN CALCULATION 
OF THE INCIDENT FLUX ON SHELL 2 

Here we show that the integrand function in equation (6) for 
the incident flux from shell 1 on point P in shell 2 comoving 
frame is insensitive to sin 2 a" so that the integrand can be 
taken out of the integral as a constant. We will also show 
that the integrand is of order unity. 

Let us call the integrand function F(a, 8). All equations 
we have are 

(1 + /3 2 cosa) 3 



F(a,l 



tana 



\J3i + cos(a - 61)] [1 + fa cos(a - 6>)] 2 



r (1 Pa < Post's 

12(1 - -3- + -5- J 



tan a = 



r 2 (cosa + /3 2 )' 



(A-l) 
(A-2) 

(A-3) 



We want to precisely estimate F(a, 8) and its dependence on 
sin 2 a" in the range of a" from to a", where the subscript 
"j" always denotes the edge of shell 1. Through equations 
(A-2) and (A-3), we can express 8 in terms of a only: 



sin a(l 



§2 , fe«N 
fl, ~T R, J 



cos a + fa 



(A-4) 



Then F(a,8) becomes F(a). Also, from equation (A-3), ex- 
press sin 2 a" in terms of a only: 
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sin 2 a + (cos a + fa ) 



(A-5) 



Therefore, we can plot F(pt) numerically as a function of 
sin 2 a". 

The only thing left is to calculate the upper limit of a. 
Equation (A-2) gives 



tana./ 



1 2 (1-7J7 + — 



(A-6) 



Denote k = tana", and apply it onto equation (A-3) and 
square both sides of (A-3). Then we get a quadratic equation 
of cos Oii 



1 2 P 2\ 



(1 + fc 2 r 2 ) cos 2 ay + 2/3 2 fc r 2 cos ay + A; F 2 /3 2 -1 = 0, (A- 7) 
with roots 



cos a, = 



-fak 2 Tj± yw+T 
1 + fc 2 r 2 



(A-8) 



The second root (— ) can be ruled out, because when we put 
it back into equation (A-3), the second root gives k < 0, 
while equation (A-2) requires k > 0. Thus 



a,- = arccos 



-faV\k 2 + yTO 

1 + fc 2 r 2 



(A-9) 



is the sole root of the upper limit of a. 

We numerically plot F(a) vs. sin 2 a" for the following 



model parameter space: ^ ranges from 0.05 to 10, - 

ranges from to 1000, and 9j= 0.1. We find, for ^ < 1, 
.F(a) is approximately a linear, monotonically decreasing 
function of sin 2 a". Its maximum is 1 and is at a = 0, and 
minimum is always > 0.01 and is at a = aj. A smaller £2 
or a larger St always gives a smaller aj , thus a minimum of 
F(a) closer to 1. 

This clearly shows that for -j^ < 1, F(a) is very weakly 
dependent on sin 2 a" thus can be taken out of the integral 
as a constant ~ 1. 

For ^ > 1, F(a) is not always on order of unity. Since 
shell 2 is moving toward shell 1, the relativistic beaming 
effect is important (Panaitescu 2007) . When 5t is very small, 
shell 2 is very close to shell 1, so that a" is larger than j^— , 
where F re ; — is the relative LF between two shells. In 
this case we find that at an angle a" < a", F(a) starts 
to drop very sharply from on order of unity to infinitely 
small, and that dropping-down angle is approximately equal 
to tt~ — . This means that, in the case of E 2 - > 1, sin 2 a'l in 
equation (14) should be replaced by a more accurate term 
min(sin 2 a'', -pj— )• However this is only required when 5t is 

rcl 

small so that - 3 ,7„?i, < 100. Thus the more accurate flux 
ratio curves for > 1 in our Figure 2 will be slightly flatter 
in the region of - ' < 100 than the ones shown. But 
our conclusion about a faster shell 2 is not affected. 
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